Tunneling in double barrier junctions with 'hot spots' 

D. Herranz, F. G. Aliev 
Departamento Fwica Materia Condensada, 
Universidad Autdnoma de Madrid, 28049, Madrid, Spain 

C. Tiusan, M. Hehn 

Institut Jean Lamour, Nancy Universite, 
545O6 Vandoeuvre-les- Nancy Cedex, France 

V. K. Dugaev 

Department of Physics, Rzeszow University of Technology, 35-959 Rzeszow, Poland 

J. Barnas 

Department of Physics, Adam Mickiewicz University 6I-6I4 Poznan; 
and Institute of Molecular Physics, 60-179, Poznan, Poland 

Abstract 

We investigate electronic transport in epitaxial Fe(100)/MgO/Fe/MgO/Fe double magnetic tun- 
nel junctions with soft barrier breakdown (hot spots). Specificity of these junctions are continious 
middle layer and Nitrogen doping of the MgO barriers which provides soft breakdown at biases 
about 0.5V. In the junctions with hot spots we observe quasi-periodic changes in the resistance 
as a function of bias voltage which point out formation of quantum well states in the middle Fe 
continuous free layer. The room-temperature oscillations have been observed in both parallel and 
antiparallel magnetic configurations and for both bias polarizations. A simple model of tunneling 
through hot spots in the double barrier magnetic junction is proposed to explain qualitatively this 
effect. 
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Recent theoretical predictions followed by experimental observations of coherent 

tunneling in magnetic tunnel junctions (MTJs) with MgO barriers js-?] have busted re- 
search of F1/I/F2 MTJs 8| (here Fi and F2 are ferromagnetic layers and / - insulating 
barrier), and opened new perspectives of their applications in spintronic devices. For thin 
enough electrodes, electron tunneling may reveal resonant features due to quantum well 
states (QWS). Early studies explored the simplest way to realize resonant tunnehng by 
growing a thin nonmagnetic layer (NAd) between ferromagnetic electrode and barrier in stan- 
dard MTJ structures F1/NM/I/F2 |9l-lll|. However, resonant tunneling in double MTJs 
(DMTJs) Fi/ 1 / Fc/ 1 / F2 {Fc is central layer) may have advantages in comparison with the 
standard MTJs, mainly due to their enhanced tunnelling magnetoresistance (TMR) 



12 



14| and resonant spin-torque effects |l5l-ll7l|. TMR in DMTJs only weakly varies at low 



bias voltages [18|, which is crucial for applications. Last but not least, the current driven 
magnetization reversal in DMTJs occurs at relatively low current densities 17|. 

In DMTJs the QWS can strongly influence electron transport only if the Fc layer has 
submonolayer roughness and its thickness exceeds Inm minimizing Coulomb blockade (CB) 



effects 



19| . On the other hand, the Fc layer should also be thin enough so that energy 



separation of QWS substantially exceeds the thermal energy. These conditions are hardly 
fulfilled in the macroscopic DMTJs 



20 



2l| . where resonant tunneling was not observed 



mainly due to the absence of the atomically flat surfaces over entire junction lateral dimen- 
sions. Evidence for a local tunneling through QWS in the central Fe layer mFi/I/FjI/NM 
junctions was provided by lovan et al. 22| using the point contact technique. Recently, 



Nozaki et al. 



23| reported on resonant tunneling effects in macroscopic DMTJs with Fe 



nano islands incorporated into the thick MgO barrier. Their results have been interpreted 



as due to combined QWS 



2^ and CB effects. 



This letter reports on the detailed study of electron transport in epitaxial macroscopic 
Fe(100)/MgO/Fe/MgO/Fe junctions with continuous middle Fe layer and current flowing 
through the 'hot spots'. Oscillatory conductance and TMR with applied voltage present 
clear signatures of the local coherent tunneling through QWS. The oscillations have been 
observed in both parallel and antiparallel magnetic configurations and for both bias polariza- 
tions, contrary to Ref. 23| where they were seen only in the parallel configuration and for 
one current orientation. Moreover, we observe oscillations in the room temperature (RT) 



regime, while in Ref. 



23| they were seen mainly at low temperatures. Finally, the technique 
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we used for their observation is also different and makes use of some features of breakdown 
junctions, which allow to observe quantum oscillations in a continuous layer. Our results 
have been explained within a simple model which assumes formation of single or multiple 
'hot spots'. 

The junctions under study have the following structure: MgO/ /MgOionm /Vi.5„m /Cr^onm 

/Co^nm /Fesnm /MgOznm /Fe^nm /MgOsnm /FCiOnm /COsOnm /PdlOnm /AuiOnm- The MTJ 

stacks were grown by molecular-beam epitaxy with the base pressure of 5x10^^*^ Torr in the 
presence of atomic Nitrogen. The MgO barrier and the Fe layers were grown at RT. The Fe 
was annealed to 450° for flattening. The high resolution cross-sectional TEM images (Fig. 
la) generally corroborate the good structural quality and homogeneity of DMTJs. Since the 
TEM images are deduced from inverse Fourier transform of a diffraction pattern integrated 
across the sample thickness, the amorphous zones, if present, add a diffuse backgound to 
the diffraction pattern of the single crystal zones. More details of standard double barrier 
samples growth may be found in j20|. The specificity of the samples studied here is a 
Nitrogen doping of all the layers achieved during the growth, with a Nitrogen concentration 
roughly estimated to less than 2%. Although further studies are needed to determine the 
concentration of Nitrogen inside MgO, recent report 25| revealed that MgO barrier in 
MTJs may be doped up to 2.5% of Nitrogen without changes in the crystalline structure. 
In our samples, the structural analysis by RHEED, Auger spectroscopy and magnetometry 
demonstrate that, despite the Nitrogen doping (as evidenced from Auger, see Fig. lb), the 
structural and magnetic properties are not affected. The RHEED patterns of Fe and MgO 
(not presented here) are identical to those of Nitrogen-free samples. What is important 
for the studies presented in this paper is that, from electrical transport point of view. 



the barrier doping by Nitrogen is responsible for local 'soft' dielectric breakdown [26| with 
reduced breakdown voltage. This may be expected to keep barrier and central Fe electrode 
compositions nearly unchanged. After the MBE growth of the multilayer stack, the MTJ 
structures were patterned to lOxlO/xm^ by UV lithography and Ar-ion etching, step-by-step 
controlled in situ by Auger. 

Figure 2a shows the RT magnetization measured with magnetic field along the Fe easy 
(100) axis (EA). One observes well defined transitions of the three distinct Fe layers. This is 
indirect indication of epitaxy and conservation of magnetic properties (fourfold anisotropy) 
for Nitrogen doped Fe. The zero bias TMR is close to 30% (Fig. 2b). 
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Application of the bias exceeding roughly 500 mV leads to breakdown of the junctions, 
which in turn reduces zero-bias TMR down to about 4% (see Figs. 2c, d) . Two other obser- 
vations indicate indirectly that this breakdown decreases the effective MgO barrier height 
most probably by transforming locally nearly crystalline MgO regions near a 'hot-spot' into 
amorphous. Firstly, the coercive field of the central free Fe layer remains unchanged after 
breakdown (Hc=45 Oe), while coercive fields of the upper and bottom electrodes increase 
substantially about 50%. These changes are most probably due to hardening of the Fe/Co 
interfaces by high current density close to the 'hot-spot'. The second indication for the pos- 
sible amorphization of the MgO barrier during the breakdown with intact middle electrode 
is our experimental observation of the signature of QWS in the electron transport in some 
of the broken DMTJs which we discuss below. 

We have carried out detailed study of the RT resistance as a function of magnetic field 
with three different field orientations: EA, the hard (110) axis -HA, and in intermediate axis 
(I A) situated approximately at 10 degrees from the EA. Here we present results obtained 
at bias voltages up to 1.5V and with the steps of 25mV for those broken DMTJs which 
showed reproducible signatures of the changes in the resistance with bias (Fig.3). The zero 
bias TMR varies substantially along the three mentioned directions. As expected, TMR 
is reduced along lA in comparison with EA and HA directions, where transitions in all 
three magnetic layers as indicated by arrows are reflected in the TMR. We shall concentrate 
further on the electron transport data obtained with field along HA and lA directions, 
where the largest relative changes in TMR vs. bias were observed. 

Figure 3a shows typical bias dependence of the resistance R for parallel (P) and antipar- 
allel (AP) states with magnetic field along HA. One observes oscillatory behavior of R and 
TMR with a period close to 150mV in both P and AP states. It is important to note that 
these oscillations, more clearly resolved for negative bias when current flows from the upper 
to bottom electrodes (Figs. 3a, c), have period that is in reasonable agreement with the 



predictions by Wang et al. 



24l |. Figures 3a, c mark with arrows the majority spin QWS 



energies calculated within about IV above (red) and below (black) the Fermi energy 24l |. 

Although the absolute values of TMR measured with the field along lA are reduced in 
comparison with those for the fields along EA and HA, the relative changes of the TMR with 
bias are substantially enhanced (see Fig. 3b,d)). In order to understand this effect we reminc. 
that the measurements of tunnel resistance with the field along I A are usually observed {27 1 
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to be most sensitive to small variations in the angle between magnetizations of the fixed 
and free layers in comparison with EA and HA configurations. We suggest here that the 
strongest relative changes in TMR(V) for the lA configuration could be a consequence of 
local spin-torque effects which are predicted to be enhanced with intermediate alignment of 
the ferromagnetic layers [3]. Figures 4 a, 6 represent 3D plots of TMR vs magnetic field 
and bias with magnetic fields applied along HA and lA respectively. Dependence of TMR 
on bias is observed to be more asymmetric with field along I A (Figs. 4 a,b), which is in 
agreement with possible infiuence of local spin torque effects in the breakdown regions. 

Let us now discuss physical mechanisms which could be behind the main experimental 
findings. Before breakdown, the current is roughly uniform across the junction area and 
weak interface disorder might introduce decoherence suppressing the effects due to QWS. 
The 'soft' breakdown of the doped MgO may create defects and local amorphization which 
locally reduce the MgO barriers leading to 'hot spots' which connect Fe and Fe/Co leads 
with the central Fe layer. With 'hot spots' of sufficiently small lateral dimensions, electrons 
tunneling to the central layer can sample well defined structure of QWS due to lack of 
decoherence. The remaining part of the macroscopic DMTJ provides then some averaged 
featureless background signal. The feasibility of the above scenario is qualitatively supported 
by the good correspondence of the observed periodic variations in R vs.V and the theory 
taking into account QWS formed within 4.6 nm thick central Fe electrode [24] (see Figure 
3a,c). 

In order to describe the observed features, we consider a DMTJ (Fig. 4c, d) with barriers 



including a number of 'hot spots'. The average conductance of the structure is a 

J Pi(Ti) ai(Ti) dTi,wheYe ai(Ti) is the conductance due to a single i-th spot, Pi(Ti) is 
the probability of realization of a certain configuration of the i-th spot, and Fj is a set of 
parameters characterizing this configuration. Let us assume that a particular 'hot spot' is 
characterized by its lateral dimension a. The current through the spot can be then calculated 
as 

y (27r)3 % 2m 

X ^ ^kik2('^) ^fcife2^fc2fc3^35 (1) 

k2<k2m 

where fci, ki^ are the normal and lateral wavevector components of the incoming wave (layer 
1), /c2 and are the normal components of the wavevectors in layers 2 and 3, respectively. 



Here we assumed that the in- plane component of k is conserved for tunnehng from layer 2 to 
3, whereas for tunneling through the spot there is no conservation of the in-plane component 
due to broken translational symmetry, and the scattering can be described by an angle 
distribution function WkikalQ-)- In our calculations wc use the approximation WkikzlQ-) — 
^2g-afc2!^ which means that scattering to the state with large in-plane component k2i is 
effectively suppressed. Equation (1) includes transmission probabilities Tk^^Q and T^^ks for 
tunneling from the layer 1 to 2, and from 2 to 3, respectively. The sum over k2 runs over 
discrete values satisfying the quantization condition — mr. It should be emphasized 
that this condition is related only to thickness L of the layer 2 and is the same for any other 
'hot spot'. 

Calculating the integral over ku we find 



The conductance I /V as a function of bias V for a single spot is presented in Fig. 4e for 
different values of a. We assumed L = 4.6 nm, the barrier width Lb = 2.4 nm, ep = 0.9 eV, 
and the barrier height Ub = £f + 3.8 eV. As we see, the oscillation peaks related to the 
level quantization in the layer 2 are more pronounced for wide spots, and they are effectively 
damped for small a. This is because the small spot enables tunneling with nonconserved 
in-plane component of the wave vector. 

Taking into account tunneling from many different spots, we obtain qualitatively the same 
picture corresponding to a mean value of a, and proportional to the number of spots. The 
total conductance of the structure includes a constant non-oscillating part, cxo, related to the 
tunneling without spots. In Fig. 4f we present resistance calculated as R— {ao + Niai{a))~^, 
where Ni is the number of spots and a is the mean value of a. We note that the variation 
of barrier heights at the spot does not affect the position of peaks and does not change the 
shape of peaks, changing only the amplitude. Thus, averaging over randomly distributed a 
and barrier heights Ub gives a picture like for a single spot with some mean values of a and 



The applicability of the above model requires a number of conditions. Although location 
of the oscillation peaks does not depend on the spot dimension (a), the 'hot spot' should 
not exceed average dimensions over which central electrode is atomically flat. Also, as seen 
from Fig.4e,f, decreasing the parameter a makes oscillations less visible (damped). Therefore 





Ub. 
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we expect the hot spots to be roughly of nm lateral size. Secondly, if magnetic moment 
of one or two layers is reversed (AP alignment,) the resistance becomes very large because 
the minority Ai band will be displaced well above the Fermi energy. However realistically 
speaking the hot spot region of the MgO barrier may be not fully epitaxial either due to 
structural defects and/or Nitrogen doping. This can strongly reduce spin filtering by the Ai 
band in the AP alignment, substantially suppressing TMR in real DMTJs with 'hot spots'. 

In conclusion, we have presented evidence for local resonant tunneling through quan- 
tum well states in the middle continuous free layer of a double magnetic tunnel junction. 
The oscillations have been observed at room temperature in both parallel and antiparallel 
magnetic configurations and for both bias polarizations. Owing to specific features of the 
breakdown junctions, we were able to observe quantum well states in continuous magnetic 
layers. To observe similar effects in the Nitrogen free DMTJ, the junction area should be 
smaller than the size of terraces of the Fe central layer. Understanding of electron transport 
in magnetic tunnel junctions with defects and 'hot spots' is of great importance both from 
fundamental and applied points of view. Recent reports link spin torque oscillations with 
record low bandwidth to the presence of defects and hot spots inside the MgO barrier of 



MTJs 



29|. 
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I. FIGURE CAPTIONS 
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Energy (eV) 

FIG. 1: (a) Cross-sectional TEM image of the DMTJ. (b) Auger data which shows presence of 
Nitrogen doping in the Fe and MgO layers of DMTJs 
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FIG. 2: (a) Magnetization curve for unpatterned DMTJ. Arrows indicate magnetization config- 
urations of top (black), middle (green) and bottom (red) layers, (b) Typical zero bias TMR of 
the DMTJ without hot spot measured along the EA. Figures (c) and (d) show correspondingly 
zero-bias TMR for field along the EA, and the HA in junctions with hot spots. The vertical 
arrows remark coercive fields. 
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FIG. 3: Resistance vs. bias for P (black line) and AP (red line) states with magnetic field along 
the HA (a) and lA (b). Here we assign as AP state the one just after inversion of the central Fe 



24| 



layer (green arrow in Fig. 2). The red and black arrows indicate predictions by Wang et al. 
for the resonant tunneling in the parallel state with QWS above (red) and below (black) the Fermi 
level. Parts (c,d) show TMR vs. bias for magnetic field along the HA and lA. 
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FIG. 4: 3-D-plot with magnetic field along x, bias voltage along y and TMR along z directions. 
Part (a) corresponds to magnetic field along the HA and part (b) along lA. (c) Schematic pre- 
sentation of the model with a single spot and (d) the corresponding energy profile, (e) Calculated 
conductance {I /V) for DMT J with a single spot in P state for different spot dimensions, (f) R vs. 
V for DMT J with multiple spots in the P state for different average spot dimensions. 
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